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ABSTRACT 

The existing analytical solutions for gas lubricated stepped 
sector thrust bearings are reviewed. An experiment performed at high 
bearing numbers (123 to 165) using air at atmospheric pressure is 
described. Results are presented of tests conducted using a bearing 
optimized by the methods of Ausman and Wildmann for high speed operation. 
Most of the design optimizations are confirmed for high bearing numbers; 
the exception is load carrying capacity. The load capacity curve ob- 
tained takes the same form as that predicted by theory but is lower for 


, - ¥ 
any given bearing number by a constant percentage. This points to a 


need for re-evaluation of determination of load carrying capacity. 
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NOMENCLATURE 
Bearing area Gem) 
Radius ratio = r,/R 
Lubricating film thickness (in.) 
Film thickness in leading edge region (in.) 
Film thickness in trailing edge region (in.) 
Load parameter = W/p, A 
Number of sectors 
Rotation rate (rpm) 
Pressure (psi) 
Compression ratio = h,/h, = (A a DH = 1+ 1l/~ 
Mass flow rate (lbm./sec.-in.) 
Radial coordinate (in.) 
Bearing inner radius (in.) 
Bearing outer radius (in.) 
Time coordinate (sec.) 
Velocity of fluid (in./sec.) 
Bearing load (lbf.) 
Axial coordinate (in.) 
Minimum clearance ratio = ho/A 
Step depth (in.) 


Angular width of groove between thrust pads 
(degrees or rad.) 


G 4) R? 
Pa he 


Bearing number based on step depth = 


Bearing number = 


64.0 R* 
pear 


Modified bearing number =f. (1-) 


Modified bearing number based on step depth = N'(i-b*) 


Sector arc length (degrees or radians) 
Step ratio = t/y 

Lubricant density (lena/ane) 
Coordinate angle (radians) 


Viscosity (reyns = pena ane 


Xr -— O- sa < 


Arc length of leading edge region (degrees or radians) 


w Angular velocity (rad./sec.-) 
Subscripts 

i Refers to leading edge region 
2 Refers to trailing edge region 
a Refers to ambient condition 

r Radial 

8 Circumferential 

A Refers to the step 

opt Optimized 
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i Introduction. 

In a rapidly expanding technology, which creates demands for machin- 
ery components that will operate in progressively more severe or unusual 
environments, it is soon found that "conventional" lubricants for moving 
or rotating parts will no longer serve. Because they possess several 
unique advantages, gases are a logical consideration as a lubricant when 
seeking solutions to many design problems. 

Gas lubricants are particularly attractive at extremely high or at 
very low temperatures. Since, in contrast to liquids, the viscosity of 
gases is relatively constant and even increases slightly with temperature, 
gas lubricants have an obvious application in very high temperature rotat- 
ing machinery. At the other end of the scale, elements such as nitrogen 
and helium retain their desirable lubricating qualities at cryogenic 
temperatures. Indeed, the ''boiloff"' from liquified gases is often used 
as the lubricant for expansion turbine gas liquifiers. 

Another important advantage of gas lubricated bearings is the extreme- 
ly low power loss from friction at very high speeds. This, plus the fact 
that this power loss is relatively unvarying, dictates the use of gas 
bearings in precision applications such as high speed gyroscopes. The low 
friction loss at high speeds plus high temperature operation makes gas 
lubricated bearings the most desirable type for gas turbines and turbo- 
jet engines. 

In many areas gas lubricants can solve the problem of contamination. 
When pumping, compressing, or expanding high purity gases or refrigerants, 
the use of the process fluid as a lubricant will prevent mixing of the 
fluid with oils or greases. In food processing machinery, elimination of 


petroleum lubricants prevents contamination of foods. In nuclear reactors, 
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it is desirable to eliminate oils, greases, and solid lubricants which 
break down or form radioactive isotopes when subjected to radiation. 

As with conventional bearings, gas lubricated bearings fall into two 
major classes, self-acting, or hydrodynamic bearings, and pressurized, 
or hydrostatic bearings. There exists a wide variety of both journal and 
thrust type bearings in each class. Self-acting bearings are, in general, 
the least complex, requiring neither separate gas supplies and piping nor 
fluid passages and orifices to complicate the design of a machine. 

Unfortunately, gas lubricated bearings have only recently emerged 
from the realm of laboratory curiosities. Currently most types are sub- 
ject to serious instabilities and require extremely close machining 
tolerances with very smooth finishes. Air film thicknesses are so small 
that thermal distortions are a major cause of failures. To gain a better 
understanding of the capabilities and problems of gas lubricated bearings, 
this project was undertaken. 

Because of an interest developed while the author was working at 
AiResearch Manufacturing Company on a summer industrial tour, this study 
will be confined to self-acting thrust bearings of the stepped sector 
type. There are several arguments which support this choice. 

Because the stepped sector thrust bearing (See Figure 1 for typical 
configuration) is a common, popular type noted for its simplicity and 
relative ease of manufacture, a large volume of definitive theory and 
design data has been developed. The most recent and acceptable theory as 
yet lacks complete experimental verification. The recent development, 
at AiResearch, of a thrust bearing test rig capable of very high speeds 
affords an opportunity to obtain data in regions of operation previously 


unattainable. Arrangements were made for the construction of a second 
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Thrust Bearing Test Rig ML1606 at the Naval Postgraduate School in 


order to conduct this research. 
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2s Theoretical Considerations. 

In 1886 the English scientist, Osborne Reynolds, published the de- 
fining differential equation for thick film lubrication. Several years 
later, Kingsbury in the United States and Michell in Australia demonstrat- 
ed the applicability of the Reynolds Equation to thrust bearings by 
development of the tilting pad bearing. 12 1918 Lord Rayleigh [1] found 
that for an infinitely wide, incompressible film, the optimum shape for 
a thrust bearing is the stepped configuration. However the step bearing 
can be optimized for only one expected operating speed or load but the 
tilting pad bearing automatically adjusts its pad angle to a new optimum 
as speed or load changes. Additionally, when side leakage is taken into 
account, the stepped bearing's load capacity markedly diminishes and its 
advantage over the tilting pad bearing disappears. For these reasons 
this bearing type was little used for liquid film lubrication. It was 
not until 1950 that Archibald [2] published the first analysis of a 
finite stepped thrust bearing using an incompressible film. 

When gas is used as the lubricant, the step bearing becomes desir- 
able for two reasons. First, for a tilting pad bearing to retain its 
continuous optimum operation, the position of the pad pivot point must be 
continually changed. This is the effect of compressibility on the bear- 
ing design. A pivot is required that consists of a cam surface that is 
very difficult to machine. Secondly ‘very close machining tolerances and 
finishes are required on gas bearings and since the tilting pad bearing 
is composed of a number of close-fitting parts, it is extremely costly 
to produce. In comparison the step bearing is relatively cheap. 

When the analysis takes into account compressible lubricants, the 
differential equation becomes nonlinear and several assumptions become 


necessary in order to obtain an analytic solution. Even though the 
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resulting solution is approximate it does establish a base for optimiza- 
tion of bearing geometry and permits correlation of theory to specific 
test bearings. 

Figure 1 depicts the geometry of a typical stepped sector thrust 
pad. The upper thrust plate (plane surface) rotates with constant angu- 
lar velocity © while the lower (stepped) surface is stationary. Radial 
grooves of angular width 6 admit gas at ambient pressure to the lead- 
ing and trailing edges of each pad. The pad itself, of sector arc 
length Y¥ , is divided into two sectors, one sector of length % being 
recessed by an amount 4 relative to the other. The larger clearance 
in the recessed region is designated h, , while the smaller, or mini- 
mum, clearance is designated h, . 

In cylindrical coordinates, the compressible lubrication equations 


for slider bearings take the form [3]: 





—otsz= wy (1) 
et * 6.88 
~~ ry. 

A - 2 
Gi Or ( dz (2) 


aM 4 2 (ev,)+455 (eVe)= 0 @) 








Q, , dQ-, 1 3Q 
+ “+ — £=0 (3a) 
p dr ar. a9 
ae (4) 
ne 
where i and V, are the radial and tangential components of lubri- 


cant velocity, Q, is the radial mass flow rate per unit circumferential 


width and > is the circumferential mass flow rate per unit radial width. 
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To obtain these equations, the following assumptions were made: 
(a) the flow is steady (does not vary with time); (b) the film thickness 
is at least two orders of magnitude smaller than r or r@ =; (c) the 
axial, z, component of velocity and its derivatives are negligibly small; 
and (d) the gas lubricant behaves as a perfect gas and compresses and ex- 
pands isothermally within the bearing. 


By introducing the boundary conditions: 


V.(z=0) = 0 
V.(z=h) = 0 
V,(z=0) = rv 
V,(z=h) = 0 


and integrating equations (1) and (2) the following expressions for the 


mass flow rates can be obtained: 








3 
oe _e@h dp (5) 
12m dr 
Q,= -|- - SE (6) 
“ fa fat in? JO 
Ausman [4] introduced equation (4) and made the assumptions that (a) 
Pa may be substituted for © , ignoring the effect of compressibil- 


ity on the pressure gradient terms and (b) ai in equation (6) may be re- 
placed with ro? where r, is some constant such that bR <r, < R where 

b is the ratio of inner to outer bearing radii = r; /R. The first assump- 
tion linearizes the resulting differential equation while the second makes 
the variables separable. Further, introduction of two dimensionless vari- 
ables, 


eeauben. Ps 
P " 


>I 


= 
R 
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and substitution of equations (5) and (6) into equation (3a) results 


in the linearized, separable differential equation 





Jf. (ee ei « JP 
seer dette) Fe = AS al 
Guar R* 
where /\ — = (8) 


Equation (8) defines the dimensionless parameter known as bearing number, 


ma , The form of the solution is 


P = X (F)-Y(e) (9) 


and results in an infinite series solution which must be solved indepen- 
dently for each flow region. 

Boundary conditions state that the pressure equals the ambient pres- 
sure at both inside and outside bearing radii, and that the pressure and 
flow rate are equal at the step when approached from either direction. 
These latter conditions are met by an iterative process which seeks a 


value for r, which, when substituted into the two series solutions, re- 


sults in equal solutions at the step for both regions. 

The theoretical load carrying capacity for a bearing is found by 
integrating the expression for p over the bearing area. 

Ausman then optimized bearing geometry by independently varying 
step location, Tt/yY , and number of pads, n, and solving for these para- 
meters at maximum load capacity. He presented his results in the form 
of one table and one graph from which a designer can determine the com- 


plete optimum geometry for a given bearing size and expected operating 


point (bearing number, A). In his optimization he assumed that the 
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optimum compression ratio, q, was 2.0 where 4 = nae = (h,+ A)/h, 
=l1t+t Lhe ; This, he found from trial solutions, was a good assump- 
tion over a wide range of bearing numbers. Note in Figure 3 (which is 
described later) that the solution for constant minimum clearance is maxi- 
mum at a clearance ratio of 1.0 which corresponds to q = 2.0. 

Wildmann, et. al., in a comprehensive study [5], obtained exten- 
sive experimental data at low bearing numbers which essentially con- 
firmed Ausman's theory in the regions the study undertook to investigate. 
Ausman's results were significantly simplified by the introduction of a 


modified bearing number 


1 b4o RY 2 
A= Ae (Ib) = A (1-8) (10) 


This enabled him to eliminate the radius ratio, b, as an independent 
parameter and to present load capacity for optimum compression ratio 
(q = 2) as one curve plotted against bearing number. This curve is 
shown as Figure 2. 

As a result of experimental data, slight modifications were made to 
Ausman's theoretical design parameters. New design parameters were pre- 
sented in graphical form so that design of a specific bearing could be 
undertaken without the necessity of extrapolating Ausman's tabular re- 
sults. The first two design charts are similar to Figures 2 and 3 but 
reflect Wildmann's correlation of theory. The remaining three are plots 
of Optimum sector ratio, ic = Ty » and optimum sector arc length 
(degrees) and optimum number of sectors, n, all plotted vs. modified 
bearing number. The test bearing investigated in this study is designed 


using the parameters given by Wildmann. 
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One additional curve, computed from Ausman‘s theory and presented 
ye P 


in Wildmann's study, is shown as Figure 3 because of its usefulness in 


comparing experimental results with theory. This consists of a plot of 


minimum clearance ratio, h2 /A_ , vs. 


ew Load 
Le eee Oe Optimized Load for q = 2 


for a constant step depth, A 


Sy Description of Apparatus 

Trust Bearing Test Rig 

The primary research tool used in this experimental investigation 
is the AiResearch Thrust Bearing Test Rig, ML 1606, an extremely versatile 
device designed to measure the operating characteristics of gas lubricat- 
ed thrust bearings at very high rotation rates. Rotating speed and bear- 
ing load can be continuously varied and accurately measured while instru- 
mentation can be provided for measurement of air film thickness, torque, 
and temperatures at selected points. Design variations permit testing 
thrust bearings up to six inches in diameter and test bearings may be of 
any type or configuration. 

Photographs of the Thrust Bearing Test Rig and its associated instru- 
mentation are shown in Figures 4 and 5. A schematic of the arrange- 
ment of major components is presented in Figure 6. The rotating assembly 
consists of a hardened stainless steel shaft one inch in diameter with an 
integral 24% inch diameter thrust disc on the lower end and a small air 
driven turbine wheel on the upper end. This shaft is supported vertical- 
ly in two AiResearch ARPF compliant journal bearings. It is supported 
from beneath by the test bearing and separation of the back of the thrust 
plate from the upper body of the machine is maintained by an externally 
pressurized air thrust bearing. This configuration is capable of rota- 
tion rates up to 105,000 RPM after running through the half-frequency 
whirl region encountered in the range between 20,000 and 50,000 RPM. 

The test bearing is supported beneath the thrust runner on the top 
face of a heavy steel piston which moves vertically and through a small 
rotation angle within a cylinder equipped with a hydrostatic air bearing. 


In the tests performed the test bearing was recessed into the face 
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of the piston and was gimballed by placing it on a 3/32 inch thick sheet 
of gasket rubber. Restraint in torque was provided by securing a .020 
inch x % inch piece of piano wire in a groove on the piston face in such 
a way as to allow the end of the wire to extend into the bottom of one 
of the air inlet grooves on the test bearing. A photograph of the bear- 
ing secured in place is shown as Figure /7. 

Air Supply and Control System 

A schematic of the air supply system is shown in Figure 8. The 
system utilizes a 50 HP Westinghouse 220 VAC, 3 phase, 60 hertz, Induc- 
tion Motor to drive a Chicago Pneumatic, Model P2CB Two Stage Air Compres- 
sor with intercooler (max pressure = 450 psi). The output air is cooled 
by a Heliflow shell and tube heat exchanger, Type 9AA14S. Water and oil 
are removed by a Wright-Austin Co., separator and float trap. These are 
followed by a Morse Diving Equipment Co., air filter and a standpipe con- 
taining a dessicant to further dry the air. Airflow to the receivers is 
controlled by a Bailey Indicator Controller, Type K47, actuating a Bailey 
Diaphram Valve Type 4% XSDS-9M. The compressed air is stored in two 
American Pipe & Steel Corp., Vertical Air Receivers (each with a volume 
of 117 pe, » max. pressure = 500 psi). Air is drawn from the top of 
both receivers and is distributed at outlets in a small manifold near the 
test setup. 

Figure 9 depicts schematically the air supply and flow control ar- 
rangement to the test rig itself. Figures 4 and 5 are photographs of 
the actual apparatus. Air is provided to the driving turbine through a 
one inch line fitted with a regulator and globe valve for rough speed con- 
trol. Fine control is achieved by a % in. line equipped with a regulator 


and needle valve. This line is placed in parallel around the rough control 
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line. A X% in. line with a secondary filter and needle valve provides 
120 psig air to the hydrostatic thrust bearing. Air is supplied to the 
hydrostatic piston journal bearing at 20 psig by a % in. line with regula- 
tor and needle valve. A similar arrangement provides air to load the 
piston from beneath, however an atmospheric bleed is placed in this line 
to prevent pressure buildup due to low flow rates. Bourdon type pres- 
sure gages are located as shown in the schematic. Loading piston pres- 
sure is measured by a 14 inch mercury U-tube manometer. 

Electronic Instrumentation 

Figure 10 schematically describes the electronic instrumentation 
used to measure rotational speed, air film thickness and temperature. 

Rotational speed of the test rig shaft is sensed by a Pace Engine- 
ering Co. SU-2 Pickup which is a very sensitive coil with a core of 
flexible iron straps. The coil, when placed in the vicinity of a rotat- 
ing shaft of ferromagnetic material, generates a small A.C. voltage of 
varying frequency as a result of the rotating magnetic field due to re- 
sidual permanent magnetism in the shaft. The Pace Pickup was placed as 
shown in Figure 4 and its output was displayed on the y-axis of a Hewlett- 
Packard Model 130C Oscilloscope and a Hewlett-Packard Electronic Fre- 
quency Meter Model 500A. The output of a Donner Sine Wave Generator 
Model 1202 was displayed on the X-axis of the oscilloscope and its fre- 
quency was accurately measured by a Systron-Donner Digital Counter Model 
1013D. Since rotation rate is a predetermined figure for each data run, 
the rate ic preset on the sine wave generator and the speed of the rig 
is adjusted until a stationary elliptical lisajous figure is achieved 
on the cathode ray tube, indicating that the two frequencies are matched. 
This results in a speed measurement accurate to plus or minus one hertz 


or + 0.1% at 60,000 REM. 
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Air film thickness measurements were made using Bently Proximitor 
distance detectors. The calibration and accuracy determination of this 
system is presented in Appendix A. During initial runs to check out the 
equipment it was found that heating and cooling of the distance probe 
(located in the center of the test bearing) was more extensive and more 
rapid than anticipated. This change in temperature causes a resistance 
change in the sensing coil and hence changes the output voltage for a 
given distance. Further, due to low friction, rundown time (time neces- 
sary to stop the shaft) was in excess of one minute and during this time 
the temperature and therefore the probe zero reading (thrust runner and 
bearing in contact), changed markedly. To correct this, two steps were 
taken. 

Two air nozzles made of % inch copper tubing were secured in the 
turbine outlet throat to direct air jets against the turbine blades op- 
posite to the direction of rotation. To actuate these, two spring load- 
ed, quick acting valves were installed. This system, seen in Figure 4 
stopped the shaft in an average time of 15 seconds. 

Two copper-constantan thermocouples were installed, one adjacent to 
the distance probe face and touching it, and the other in an ice-water 
bath as a reference junction. The output was recorded on a Mosely Auto- 
graf Strip Chart Recorder Model 7100B. The second channel of the re- 
corder was used as an event marker to record points at which data was 
taken. This arrangement is shown in Figure 10. 

These two corrective measures enabled the author to establish the 
probe temperature difference between the probe reading and the zero read- 
ing (thrust runner and bearing in contact) for each data point. The 
significance of this data will be shown in Appendix C - Experimental 


Uncertainty Intervals. 


655, 


Test Bearings 


It was originally intended to test several thrust bearings of dif- 
ferent optimum configurations but time and funding limitations restricted 
the tests to one bearing. The bearing was chosen to correspond to the 
point in Ausman's Table I [4] at b = .20, A\ =A, = 160. The design 


parameters however, were taken from Wildmann [5]. These parameters are: 


b= "920 

A= A. (i-b*) = 153,7 
7 = %/ = 135 

nN =3 

¥ = 118° 

Tv = 16° 


The minimum film thickness was chosen to be 250. microinches so, for 
q = 2.0, the optimum step depth, A  , is 250. microinches. The out- 
side radius, R, was chosen as 1.125 inches. The final design based on 
these parameters is shown in Figure 11. Note the air inlet groove which 
is two degrees wide at the average radius. The material chosen for the 
bearing was "Graphmo" tool steel, a graphite-molybdenum steel noted for 
its wear resistance and low sliding friction. The bearing was manufactur- 
ed to specifications by B & W Precision, Inc., Glendale, California. 
Inspection of the bearing by optical flat before running revealed 
a flatness of + 5.0 microinches across the bearing face. Profilometer 
readings by the manufacturer showed a surface roughness of 5 microinches 
(rms). The step depth was measured on a Johansson Mikrokator Comparator 
and varied between 230 and 270 microinches with an average of 250. These 
readings were taken numerous times across all three steps. The surface 
of the recessed area was somewhat rougher than the raised portion but was 


not measured. 
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4. Experimental Procedure 

Data was obtained in such a manner that load vs. air film thickness 
could be plotted for a single constant speed. This was repeated for a 
number of different speeds to obtain a family of such curves. No data 
was taken in the region where half-frequency whirl was encountered to 
insure that the data was not affected by this instability. 

Test Rig Operation 

All runs were made with the hydrostatic thrust air set at 115 to 
120 psig and the hydrostatic piston journal air adjusted to 20 psig. 
All starts and stops were performed with a test bearing load of 3.65 
pounds (3.0 inches Hg). Prior to starting, the speed for the run (re- 
volutions per second) was preset on the digital counter by adjusting the 
sine wave generator to the same frequency (hertz). The turbine was 
then started by quickly opening the rough control valve enough to bring 
the downstream pressure to about 5 psig. The bearing noise was monitor- 
ed continuously for whirl or rubbing by listening with an automotive type 
stethoscope. Rotational speed was adjusted while watching the lisajous 
figure on the oscilloscope and reading the frequency meter until the 
speed was slightly greater than the preset rate and steady. The tempera- 
ture recorder chart was then started and the bearing load slowly increas- 
ed, taking care that the temperature did not increase at a rapid rate and 
cause excessive thermal distortion. As the bearing was loaded, the tur- 
bine air flow was increased to maintain speed, 

When the load for the particular data point was reached, the rig 
was allowed to run until the temperature reached a steady state or was 
increasing very slowly. Then the fine speed control flow valve was used 
to control turbine speed until the oscilloscope trace showed an exact 


frequency match. At this point, the distance probe output was read and 
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noted, the turbine air was secured, the braking jets were actuated, and 
the bearing load was reduced to 3.65 pounds. As the probe output was 
read, the event was marked on the strip chart recorder. At the instant 
that the shaft stopped and the thrust disc and bearing came in contact, 
the distance probe output was again read and the event marked. 

This procedure was repeated for each data point obtained, and each 
data point was repeated several times to insure correctness and to ob- 
tain a statistical average. The data obtained are presented in Tables 
1 through 4. 

Bearing Failure 

On the last run for data, a failure occurred and the thrust plate 
contacted the bearing face with subsequent rapid loss of air film and 
wiping of the bearing surfaces. The failure was initiated at a speed of 
72,000 RPM and a bearing load of 14.7 pounds. Just prior to failure a 
vibration noise began and increased in amplitude even though the load 
was reduced sharply. The center of both the bearing and thrust disc 
were wiped to a radius of approximately .35 inches which would seem to 
indicate failure due to thermal distortion. This is felt not to be the 
case, however, since neither temperature nor temperature rate were ex- 
cessive at the time, and the load, while high, was no higher than for 
several other runs when the air film thickness was even less. Rather 
it is felt that the speed used was a critical vibration frequency for 
some component, most likely the rubber sheet used to gimbal the test 
bearing. The appearance of the surfaces after failure merely indicates 
that slight thermal distortion was present at the time of failure. 

Data Reduction Calculations 

The data reduction calculations used to obtain the results presented 


in the next section are described in Appendix B. 
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5. Results 

The numerical results calculated from the data obtained are present- 
ed in Table 5. 

Figures 12 through 15 are plots of bearing load vs. air film thick- 
ness for each speed tested. All data points are plotted and extreme 
uncertainty limits are shown. Figure 16 is a composite plot of the same 
data using only the statistical averages for air film thickness. 

Figure 17 is a graph of dimensionless load parameter, L, vs. the 
modified bearing number, LS fy, for a compression ratio, q = 2.0. 

Perhaps the most useful form of presenting the results is seen in 
Figure 18. The ratio of actual load to optimized load at q = 2.0 is 
shown on the ordinate logarithmically while the clearance ratio, =< = hey 
is plotted on the abcissa. By the optimized load it is meant the opti- 
mized load from the theoretical curve in Figure 2 determined for the 
same bearing number as the particular data point under consideration. A 
clearance ratio of 1.0 corresponds to'q = 2.0. All of the data points 


from Figure 17 are shown in Figure 19. 
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Gr Comparison of Results with Theory 

In Figure 19, the results shown in Figure 17 are superimposed on a 
reproduction of Wildmann's Design Chart A [5], which predicts load carry- 
ing capacity at q = 2.0 for a given bearing number, oS « Note that 
the results take roughly the same form as the predicted curve but the 
load capacity is signifigantly lower. 

In Figure 20, the curve from Figure 18 is replotted on a graph which 
shows both Ausman's theoretical prediction and Wildmann's experimental 
results. Note that the results of this study have not been normalized 
to a load ratio of 1.0 as was done in the study compiled by Wildmann. 

Two factors are immediately apparent. The first is that the slope 
of the curve obtained in the current study is essentially the same as 
Wildmann's correlation curve. Secondly, the new curve predicts a load 
capacity that is a constant fraction, 2/5 or 0.4 of that predicted by 


the design curves [5] based on Ausman's work. 
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Te Discussion of Results 


It is felt that the careful calibration of the probes used to meas- 
ure air film thickness, plus the temperature monitoring employed, re- 
sulted in reliable air film thickness measurements. Further, the repeti- 
tion of data points and use of average values reduced random error. This 
opinion is borne out by the smoothness of the load vs. film thickness 
curves and the consistency of the load ratio vs. clearance ratio curve. 

There are two sources of possible major error. The first is due to 
the near impossibility of constructing the perfect stepped sector thrust 
bearing upon which the theory is based. It was noted that the variation 
in the step depth was plus and minus 20 microinches or 12% of the speci- 
fied depth. This plus the rougher finish on the recessed area undoubted- 
ly affected the flow characteristics of the lubricant. The second seri- 
ous condition is thermal distortion. Even though attempts were made to 
allow temperature to reach a steady state, distortion was surely present 
as shown by the nature of the bearing failure which occurred. Thermal 
distortion would result in a nonuniform air film thickness; the film 
would increase in thickness as the radius at which it is measured is in- 
creased. It is not felt, however, that either of these conditions are 
responsible for the wide disparity between theory and experiment with 
regard to load capacity, particularly since the data is otherwise con- 


sistent. 


a 


o. Conclusions 

From the comparison of results with theory, it is concluded that: 
(a) Wildmann's correlation of the variation of load ratio with clearance 
ratio is correct at high bearing numbers for the bearing tested; (b) 
Ausman's theory correctly predicts the shape of the load vs. bearing 
number curve in the high bearing number region, and (c) the prediction 
of the magnitude of the load capacity is in error in this region. 

Ausman, in the preparation of his design table, performed separate 
optimizations for bearing geometry but extrapolated load capacity from 
curves computed for a central point in the table. He warned his readers 
that the accuracy of this extrapolation probably decreased with distance 
from the center of the table. The results of this test, of a bearing 
purposely chosen from an extreme corner of his table, confirm that Ausman's 


intuition was correct. 
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9. Recommendations for Further Study 

Now that it is established that there is a limit to the range of 
accuracy of Ausman's design parameters, additional bearings, optimized 
by the procedure used in this study, should be tested until the range 
of accuracy is established. 

A computer aided analysis of the load capacity equations should be 
performed over the full range of bearing numbers using the modified 
bearing number introduced by Wildmann. The results should give some 
insight into the true nature of the relationship and could result in a 
revision of the optimized design parameters. 

Other types of thrust bearings, such as the spiral groove or her- 
ringbone patterns, await testing at high bearing numbers. The AiResearch 


Thrust Bearing Test Rig is ideal for this purpose. 
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TABLE 1 
Summary of Data -- 54,000 RPM 


RUN NO, LOAD h 


(pounds) (4 in.) 
14-1 6.59 363 
14-2 6.59 366 
14-3 6.58 371 
14-4 6.59 363 
14-5 6.59 366 
14-6 6.59 370 
15-1 8.62 318 
15-2 8.62 330 
15-3 8.62 337 
15-4 8.62 320 
15-5 8.62 325 
15-6 8.62 320 
16-1 10.65 287 
16-2 10.65 287 
16-3 10.65 287 
16-4 10.65 284 
16-5 10.65 289 
16-6 10.65 294 
17-1 12.67 254 
17-2 12.67 260 
17 =3 12507 256 
17-4 12.67 258 
17-5 12.67 254 
17-6 12.67 253 
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TABLE 1 (Cont) 


RUN NO, LOAD sy 
(pounds) (4 in.) 
18-1 13.69 246 
18-2 13.69 246 
18-3 13.69 244 
18-4 13.69 Zoe 
18-5 13.69 243 
18-6 13.69 250 
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RUN NO, 


l-l1 
2-1 
2-2 
2-3 
2-4 
2-5 
2-6 
3-1 
3-2 


3-3 


3-5 
3-6 
4-] 


4-2 


4-5 
4-6 


5-1 


5-4 
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TABLE 2 
Summary of Data -- 60,000 RPM 


LOAD 
(pounds) 


no data 
4.06 
4.06 
4.06 
4.06 
4.06 


4.06 


10.65 
10.65 
10.65 
10.65 


10.65 


3/7 


(4 in.) 


442 
452 
455 
466 
456 
455 
360 
363 
370 
380 
385 
394 
320 
313 
330 
333 
337 
328 
296 
294 
294 
296 


298 


TABLE 2 (Cont) 


RUN NO, —— h, 
( 4 in.) 
5-6 10.65 289 
6-1 12.67 256 
6-2 707 257 
6-3 M2 367 256 
6-4 12.67 268 
6-5 12.67 261 
6-6 12367 260 
7-1 13.69 243 
7-2 13.69 251 
7-3 13.69 253 
7-4 no data 
7-5 13,69 253 
7-6 13.69 251 
7-7 13.69 253 
8-1 14.19 240 
8 <2 14.19 250 
8-3 14.19 241 
8-4 14.19 240 
8-5 14.19 243 
8-6 14.19 240 
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TABLE 3 


Summary of Data -- 66,420 RPM 


RUN NO, LOAD h 
(pounds) 2 
C4 in.) 
9-1 6.59 376 
9-2 6.59 380 
9-3 6.59 385 
9-4 6.59 385 
9-5 6.59 377 
9-6 6.59 374 
10-1 8.62 332 
10-2 8.62 347 
10-3 8.62 337 
10-4 8.62 339 
10-5 8.62 336 
10-6 8.62 337 
11-1 10.65 299 
11-2 10.65 305 
11-3 10.65 312 
11-4 10.65 295 
11-5 10.65 295 
11-6 10.65 298 
12-1 12.67 263 
12-2 2767 274 
12-3 12.67 264 
12-4 12.67 274 
12-5 12.67 268 
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TABLE 3 (Cont) 


RUN NO. LOAD h 


(pounds) nl 
13-1 14.70 250 
13-2 14.70 247 
13-3 14.70 240 
13-4 14.70 247 
13-5 14.70 243 
13-6 14.70 246 
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TABLE 4 
Summary of Data -- 72,000 RPM 


RUN NO, LOAD h 


(pounds) ar na 

19-1 8.62 337 
19-2 8.62 353 
19-3 8.62 346 
20-1 10.65 301 
20-2 10.65 301 
20-3 10.65 301 
21-1 12 67 265 
21-2 12.67 267 
21-3 12.67 270 
22-1 13,69 258 
22-2 13.69 260 
23-1 14.70 253 
23-2 14.70 Bearing failure. 
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Bearing Load, W (Ilbf.) 
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FIGURE 13. BEARING LOAD vs. AIR FILM THICKNESS - 60,000 RPM 
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APPENDIX A 
DISTANCE PROBE CALIBRATION 
Selection of Measuring Systems 

The thickness of the air film between the test bearing and its 
thrust runner was expected to vary between 150 to 400 microinches. It 
was felt that significant data could be obtained with a measuring techni- 
que which could discriminate within plus or minus 10 microinches during 
actual operation. Based on descriptions of apparatus found in the liter- 
ature, two types of sensors should be considered --- a capacitance probe 
which would use the air film to be measured as the dielectric between 
two parallel plates, and the magnetic, or reluctance, type which senses 
the proximity of a conductor by the extent to which it disturbs a mag- 
netic field surrounding a sensing coil. 

Although it is capable of high accuracy, the capacitance type was 
eliminated because it depends on the dielectric constant of the air film. 
It was felt that changes in temperature and moisture content would un- 
predictably alter the dielectric properties of the air film. 

Bently Proximitors manufactured by the Bently-Nevada Corp.were chosen 
because of their high potential accuracy, convenient mounting arrange- 
ments, and low cost. Since they are of the magnetic reluctance type, 
their measurements depend upon the nature «” the conducting surface being 
sensed. Variations in the properties of nonconducting materials be-= 
tween the probe and the measured surface have no measurable effect on 
the accuracy obtained. Temperature of the probe itself was found to 
have an effect, however, even though silver-wound coils, advertised as 


the least temperature sensitive, were used. 
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Enhancement of Accuracy ; 

The selected probes were intended to measure a range of 0 to .010 
inches and to resolve to + .0001 inches. To attain the necessary ac- 
curacy (+ 10 microinches), the resolution had to be improved tenfold. 

According to the manufacturer, the accuracy and resolution is direct- 
ly proportional to the stability of the input voltage. Accordingly, an 
ultra-stable power supply meeting the voltage and current requirements 
was purchased from Systems Research Corporation. Measured ripple and 
noise of their Model 3566 Power Supply was found to be 50 to 80 micro- 
volts peak to peak and the output varied no more than 1.0 millivolts in 
a 24 hour test. 

To read the output voltage accurately, a digital voltmeter, Electro 
Instruments Co. Model DVX-500, capable of measuring to the nearest 0.1 


millivolts was obtained. 


Probe Calibration Methods 

Photographs of the probe calibration rig used are shown in Figures 
Al and A2. The apparatus consists of a Johannson Mikrokator Model §S 
510E-7 Comparator altered for use as a calibrator. The Mikrokator dial 
indicator seen in Figure A2 is accurate to the nearest one millionth of 
an inch and is thus capable of establishing the actual resolution of the 
distance probes. The comparator was altered by providing a rigid clamp 
fixture to hold a distance probe over the movable table, and by making 
arrangements to adjust the position of the dial indicator head by means 
of a fine adjustment mechanism salvaged from a microscope. 

The surface to which the probes were calibrated was the surface of 
a disc made indentical to the thrust runner actually used in the thrust 


bearing test rig. This was done by cutting the disc from the same billet 
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of steel and having it follow the shaft with thrust runner through the 
same machining, heat treating, grinding and lapping processes. 

Provisions were made for checking the effects of temperature on the 
distance measurements. The ambient temperature was accurately recorded 
before, during, and after each calibration. The calibration rig was al- 
lowed to "soak" for at least three hours at any given temperature to elim- 
inate any transient expansion or contraction. 

A temperature sensitivity test was performed to insure that tempera- 
ture changes in readout were not due to mere distance changes as a result 
of expansion or contraction of clamping fixtures of the calibration rig 
itself. A probe was clamped in a vertical position with its sensing sur- 
face up. The test surface disc was allowed to rest above the probe «ace 
with a separation of .015 inches maintained by a thin sheet of Pyrex 


glass. The coefficient of thermal expansion, ™ , of Pyrex glass is 


-6 in. 
in.°F. 


glass thickness ,a , will change by 


2.0x10 So for a temperature change of 1OOF, for instance, the 


Crate \ Coaled 


(2.0x107°) (.015) (100) 


i 


3.0 microinches 
This arrangement was then subjected to heated air blown from a portable 


heat gun while the output reading was watched for deviation. 


Results 

In the temperature sensitivity test, raising the temperature ap- 
proximately 50F caused a change in reading of about 100 millivolts. This 
represents a relative distance change of about 140 microinches and con- 


clusively confirmed suspicions that the probes themselves are affected by 
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temperature. 

Calibration curve data are presented in Tables Al-A3. The resulting 
calibration curves are presented in Figure A3. Note that data was taken 
only in a small linear region of the probe's effective range since larger 
displacement measurements were not anticipated. The slope of the curves 


was the same in every case; 


microinches 
millivolt 


Slope = 1,4118 

Discussion 

From Figure A3 it is immediately apparent that although temperature 
changes affect the position of the curve, only the zero, or reference, 
point is affected. The slope of the curve remains constant over the ex- 
pected range of measurements. Accordingly, the experimental procedure 
during actual runs included the determination of a zero point (bearing 
and thrust runner in contact) immediately after each run while the rig 
and probe were still at operating temperature. Measurements were then 
taken from this zero point using the constant slope of the calibration 
curves. 

To determine data scatter and accuracy, the data was originally 
plotted on 20x24 inch graph paper with 10x10 squares to the half-inch. 
The mean deviation from the straight line plots was 0.96 microinches. 
Thus by accuracy and stability enhancements, together with careful cali- 
bration, the resulting data may be considered reliable within + 1.92 micro- 
inches, or twice the mean deviation. This is better than five times the 
necessary accuracy and fifty times better than the manufacturer's intended 


accuracy. 
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TABEE Al 


DISTANCE PROBE CALIBRATION DATA 


68,6°F 
Run No. Relative Distance Probe Output 
Cin) (mv ) 
L 400 6069 
379 6035 
ooy 6000 
296 5968 
268 3945 
204 5900 
170 5875 
LZal. 5838 
66 58C0 
2 5756 
2 392 6056 
39 6031 
320 6005 
210. 5968 
216 5953 
168 5896 
fl 5857 
63 5820 
9 B//7 
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TABLE A2 


DISTANCE PROBE CALIBRATION DATA 


75aeo 
Run No. Relative Distance Probe Output 
fin. } (mv ) 
Z 398 5943 
350 Deu 7 
302 5874 
Lon 53836 
193 S17 
130 5138 
p5 5724 
52 5695 
4 393 5941 
344 5398 
310 5872 
229 5815 
176 oy] 
114 Shae 
62 5689 
0 5648 
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TABLE A3 


DISTANCE PROBE CALIBRATION DATA 


49 ,.4°F 
Run No. Relative Distance Probe Qutput 
(in.) (nv) 

> 396 6242 
366 6219 

307 6175 

254 6134 

234 6123 

Lia 6082 

219 6037 

68 6005 

0 2990 

6 395 6207 
360 6186 

310 6154 

251 6117 

207 6088 

158 6054 

123 6028 

87 6004 

De 598) 

18 29516 

0 9942 
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APPENDIX B 
DATA REDUCTION CALCULATIONS 
The results of measurements shown in the Summary of Data (Tables 
l through 4) were calculated from the raw data. Examples of these calcu- 
lations follow. 
Since the speed measurement is indicated by a frequency (hertz) and 
1 hertz = 1] revolution per second 
speed in revolutions per minute is computed by multiplying this value by 
60. For example: 
N = 900 hz. x 60 = 54,000 rpm 
The bearing load is indicated by the height of a column of mercury 
in inches. The manometer reading is converted to psi and multiplied by 
the piston area. The weight of the piston and test bearing assembly is 
subtracted from the result to obtain the actual bearing load: 
2.036 ae = 1.0 psi 
Piston Area = 10.317 in. 
Weight of Piston Assembly = 11.65 lbs. 


Therefore, for a reading of 4.0 in. Hg, 


eee esOmeetO). 317 
2.036 


- IlV63> =28.62) fbae 

Air film thickness is determined from the difference (in millivolts) 
between the voltmeter reading at a particular speed and load and the 
reading immediately thereafter taken with the bearing and thrust disc in 
contact. This difference is multiplied by the constant slope of the 


probe calibration curve (Appendix A) to obtain air film thickness. From 


Run No. 14-1, 
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h,= (5864 - 5607) (1.4118) 
= 257 x 1.4118 = 363 microinches 
Subsequent analysis of the results is performed using the dimension- 
less parameters, "a , Lb, =, ana Load Ratio (L/L 3 ). The calcu- 
lation of these parameters is described in the following paragraphs. 


é 


a ’ ‘ a 
The bearing number, | Bs 
\ <a 7 —_— ' es ' 
where constant values are 


4 =m28 tor” reyns 


x = 1.125 ee. 
__ = Oe. 
ee e25.00% 10 iis 
b = 0.2 
and . 
= / af ' [ f a 
<h- = fm ae aa = \ 


St 
nn 


p = barometric pressure (in. Hg) x .4912 in, H5. 


For Run No. 6, 


1000 (277) = 6283 rad./sec. 


p = 30.31 x .4912 = 14.89 psia 


and 
-~10 Z 
a (o7(28 x 10 ) (6283) (1.125 ) Gig ri 
oo (14.89) (25 x 1U79): . 
= 137.8 
The load parameter, L, is given by 
W 
L=- 
P.A 
where 


A = (Ro ~ 52) rx (354/360) 


(1.1257 - 225°) 7 (354/360) 


tl 


3.7 DD in. 


75 


The fraction 354/360 accounts for the fact that the air inlet grooves 


are not part of the load carrying area. Again for Run No. 6, 


i 
(14.89) (3.755) 
= 0.227 
The clearance ratio, “& = ho/A , is computed from the average 


film thickness for a given run. Using Run No. 6, 


h, (avg) = 259.7 microinches 





and 
ee 
XY = 750 1.04 
The optimized load for q = 2 is found by consulting Figure 3 for 
/ 
the particular bearing number, Te . For example if 
A, = 137.8 
L = .628 
and 
pe eee, a 
Load Ratio = 628 361 
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APPENDIX C 


EXPERIMENTAL UNCERTAINTY INTERVALS 


Reference |9/ gives the uncertainty interval for a result as 


woaere WwW 


- 2 ‘ a as 
ae \ Tey oa Vs 7 
4 


W jes, We are the uncertainty intervals of the variables 


making up the result, R is the expression for the result, 


V 
ry 


are the variables making up the result. 


The expression for bearing load is 


where 


; \f “ \ . ~ 
load pressure (in. Hg) 
piston area (in.) 


weight of piston & test bearing assembly (lbm.) 


The uncertainty interval for bearing load is 


Dividing through by W to nondimensionalize the equation results in 


For 


W,, / 
ID.65 


iH 


= « ee > 
MW. j 3 \y 


me . enue ene © 
es ies * 5 as ae? 
ee i \ F F 2 wis 


a 


= a f i 
af “ ? We i Wan e / Why é (2 
saa, now t- r } PrP 
~ y «=~ +7 aT + ‘ a 
VAN ran / ES a \ WA a 


i” 


4.4 + .05 in. Hg (26560 1) 


10.627 @ 2081 in. (20 to 1) 


11.65 + .@2evim (20®to~1) 

aN =e, Coy fem /, a = 
“| (Ga d | a sa "Ailes 
.01544 
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and v 


> 


= 2 


and 


w= (10.65) (.01544) = + 0.164 bf. 


w7 


Air film thickness is computed in the following way: 


hy (Vi - V,)S 


where 


= 
i] 


probe output voltage at speed and load (mv) 


<= 
i] 


probe output voltage at zero point (mv) 


op] 
il 


slope of probe calibration curve CA in. /mv) 

One additional factor must be taken into account. Inspection of the probe 
calibration curves shows an average temperature sensitivity of 14 mv/°F 

at a constant distance between probe and measured surface. Examination 

of all temperature data taken shows that in only one out of ten readings 
did the probe temperature differ at speed and load from the Zero point 

by more than 1.25 °F. An actual temperature trace from the strip chart 


recorder is shown as Figure C-l. Defining an additional error, 


W 


= uncertainty interval due to a difference in prcbe 


temperatures 


(14.0) (1.25) 


ery 
It is felt by the writer that repetition of each data point an average of 
six times halved this uncertainty interval at the very least. The inter- 
val used is 

wo 8.75 mv (10 to 1) 


Therefore for Run No. 5-2 using the previous method, 


/ 
Wie " [Wey ez Ws = Wy = (| = 
he \¢ Mi - ( | i ae - ~ 





whe re 
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V, = 5788 + 2 mv (2050 1) 

V~ | 5580 + 2 mv (20.20.42) 

S = 1.4118 + 005 yin. /my (20 G) 1) 

and | V 

va | 2 —  \ eer fe )3 YP | - 
= soe ee eh _ = wm 
er ( " Sxo] \ 28 VhAHa / 

WwW = 294(.0422) 


=°12.41 microinches 


The expression for bearing number is 


f ( : \ We | ( | *) 
Oe a {- |, 
& \ 


f 
‘ ‘ (is A, 


and by similar means 


ae" : a Ys 
sf \ + Oe Ne /. - 
ee — Ca Ai “ ) + ey . a (a) | 
A,’ t ( te ‘ bo A 


4 





where for Run No. 5-2, 


A = 28x10719 4+ ipsam reyns (20 to 1) 
) = 6283 + 6.283 rad./sec. (20 to 1) 
Bo > 20651 + .01 in. He (20 to 1) 
A = 250 + 10 microinches (20 bo, 1) 
and 
- | sae 
ch! | ~ (a ‘ae _ 
1 ye [eye > Gey] 
Ww, = 137.8(.0669) 
= + 9,22 


Similarly for Run No. 5 (averaged points), 


vee, 


eg 


t 
1 








Wa | (4. ) + Wa \~ 
c<. | He A 


Woe ; all 10 yy /2. 
ie i Fe a es 


$1, 1S g@s058 h) 











= 
| 


oX, 


+ .069 

The foregoing are examples of the determination of uncertainty 
intervals for each data point. These intervals are represented on all 
eraphical presentations of data as short lines through the data points in 


the following manner: 


UNCERTAINTY 
Inte gvAcs 
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FIGURE Cl. PROBE TEMPERATURE AS A FUNCTION OF TIME 
FOR RUN NO, 17-4 
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